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Термоупругодинамическая неустойчивость решения контактной задачи для покрытия с 
учетом тепловыделения от трения* 


В. Б. Зеленцов, Б. И. Митрин, С. С. Волков, А. С. Васильев 


Рассматривается одномерная задача термоупругости о вертикальном внедрении жёсткой ПоОЛУупПЛлоскости, движу- 
Щщейся Г оризонтально с пост ОЯННОЙ скоростью по поверхности упругого покрытия (1 ПОЛОСЫ В плане), нижняя сторона 
которого жестко оперта на недеформируемое основание. На поверхности основания поддерживается постоянная 
температура. Тепловой ПОТОК, порожденный Фрикционным контактом, направлен в покрытие. Решение задачи по- 
лучено с помощью интегрального преобразования Лапласа и представлено в виде контурных интегралов. Иссле- 
довано расположение полюсов подынтег! ральной функции решения в КОМПЛЕКСНОЙ ПЛОСКОСТИ при различных зна- 
чениях параметров задачи. Распределения температуры, смещений и напряж ений по толщине покрытия получены 
в виде бесконечных рядов по собственным функциям. Показана неуст. ОЙЧИВОСТЬ полученных решений задачи на 
всём интервале времени при Любой скорости скольжения полуплоскости по поверхности покрытия. 

Ключевые слова: термоупругодинамическая неус ТОЙчИвОС 7ть, связанная задача термоупругости, Фрикционный 
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Введение. Теоретическому и экспериментальному исследованию термоупругодинамической неустой- 
чивости скользящего фрикционного контакта уделяется достаточно большое внимание со стороны 
научно-технического сообщества [1-9]. При теоретическом изучении задач динамики термоупругого 
скользящего контакта наиболее часто использующимися методами исследования являются методы ма- 
лых возмущений [3—5], с помощью которых устанавливается термоупругодинамическая устойчивость 
или неустойчивость решения задачи, определяется параметрическая область устойчивости или не- 
устойчивости решения задачи [6]. 

Постановка задачи. Для исследования динамики скользящего фрикционного термоупругого контакта 
рассматривается одномерная задача плоской деформации о скольжении с постоянной скоростью И жест- 
кой полуплоскости В (р <х< о) по поверхности (х = в) упругого покрытия в виде бесконечной по- 


лосы шириной № (@0<х <), нижняя сторона которого (х =0) жестко соединена с недеформируемым 
основанием в виде полуплоскости А (- ©<х< 0) (рис. 1) [5]. 


| 
жёсткая полуплоскость —_—_ 


ПРИ ГРИ 





упругое покрытие 





А жёсткая полуплоскость 


Рис. 1. Схема скольжения жесткой полуплоскости по поверхности упругого 
покрытия, жестко соединенного с недеформируемым основанием 





* Работа выполнена при поддержке грантов РФФИ № 14-08-91166-ГФЕН_а, 14-07-00271-а и 14-07-00705-а. 
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Скольжение недеформируемой полуплоскости В’ по поверхности покрытия происходит с учетом 
кулоновского трения, но без учета износа покрытия. Полуплоскость В деформирует упругое покрытие, 
смещаясь вдоль оси Х по закону (А, Е) =-А(Ё), #>0, где ым(р, Е) - вертикальное смещение поверх- 


ности покрытия при х=А. В начальный момент температура покрытия нулевая: 7(х,0)=7Т =0 
(0 <х< В), где 7 (х, Е) — функция распределения температуры в покрытии. Движущаяся полуплоскость 


ОТ (В, Е) 
х 


В теплоизолирована, а поток тепла К ‚ образующийся за счет трения, направлен в упругое 


покрытие. Так как нижняя сторона покрытия лежит на недеформируемом основании (х=0) в виде 
полуплоскости А, то на этой стороне покрытия упругие смещения (/(0, Ё) =0. На нижней стороне по- 
крытия поддерживается нулевая температура 7(0,Е) =0. С учётом того, что до начального момента 
ди(х,0) _ 
9 
Сформулированная одномерная задача о скользящем контакте приводит к следующим гранич- 

ным условиям: 


ди 
времени покрытие находилось в покое, начальные условия на Ш и = нулевые: и(х,0) = 0. 


х=р:  иВЮ=-МО, КОТО, - ле, #>0; (1) 


х =0: (0, Е) =0, 7 (0,2) =0, >0, (2) 
где о(х, Е) — Напряжения сжатия в покрытии, К- коэффициент теплопроводности материала покры- 
тия, А(Е) — закон внедрения полуплоскости Вв упругое покрытие, который принимается в следу- 


ющем виде: 
0, —-®<Ё<0, 
А(ЁЕ) =А.)-1+е*, ОЗЕРЕ ЕО, (3) 
1, ЕЕ <Ё<о, 
где =>0, К =='И2, Л, - глубина внедрения полуплоскости В в упругое покрытие 
(0<Л, <А). 
Скорость внедрения полуплоскости В определяется как производная А(Е) и выражается формулой 
0, —®<Ё<0, 
А(Е)=А =, — 0<Ё<Ё, (4) 
0, Е <Ё< о. 


Максимальная скорость внедрения достигается при Ё = & и составляет Д(Е.) = 2Л,= , при этом 
начальная скорость внедрения равна Д(0) = и, = До, откуда = =И,/Ду (при заданных и, и Дь). 

Изменение (/(х,Ё) и о(х,Ё) в полосе описывается одномерным уравнением теории 
упругости [10] 


2 
В О: 0<х<й, Ё>О, (5) 
ОХ ОЁ 
а температура 7(х, Е) - уравнением теплопроводности [11] 
2 
Е 0<х<й, Ё>О, (6) 
9х? коЁ 


где р, к - плотность и коэффициент температуропроводности материала покрытия соответственно. 
Связь между о(х, Е), и(х,Ё) и Т(х,Е) устанавливается соотношением Дюамеля - Неймана [12] 
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21-м) ди 2ы+у) 
1-2, ох 1-2\ 
где |, м, а - модуль сдвига, коэффициент Пуассона, коэффициент линейного расширения материала 

покрытия соответственно. 
Для удобства построения решения задачи (1) - (6) выражение ов(х, Е) (7) подставляется в диф- 
ференциальное уравнение движения упругой среды (5). В результате этого получим дифференциаль- 
ное уравнение движения упругой среды относительно ш(х,Ё) и 7(х, Е) 


2 и — 
д`и о ОА: Раб. а 2ы( м (8) 
9х? а? 0? 1-м ах Ра - 2) 
где а - скорость упругой волны. 


Начальные условия на (/(х,ё), 7(х,Ё) и А(Ё) нулевые 
ди(х,0) _ 
ОЕ 





с(х,Е) = а7(х,Е), (7) 





и(х,0) = 0, 7Т(х,0)=0, 0<х<йА, 4(0)=0. (9) 


Таким образом, решение рассматриваемой задачи, с учетом начальных условий (9), сводится к решению 
дифференциальных уравнений (6), (8) с граничными условиями (1), (2), а о(х,Ё) в (1) определяется 
формулой (7). 

Решение поставленной задачи. Поставленная нестационарная начально-краевая контактная задача 
(1)-(9) относится к разряду собственно-связанных задач термоупругости, так как о(х,Е) и 7(х,Ё) свя- 
заны не только в формулах (5) и (7), но и во втором граничном условии (1). Её решение может быть 
построено различными методами математической физики [11]. С помощью интегрального преобразо- 
вания Лапласа [13] решения поставленной задачи (1)-(9) 7Т(х,Е), и(х,Ё), о(х,Ё) определяются с 
помощью контурных квадратур 


1-иУуй 1 я с Е 
рота ЕО К(едекрЕР 2 , Г _ (10) 
и(х,Ё) = т: [22 м, (х, 2) (еек р (11) 
Г 


о(х,Ё) = я р Е, [5 (х, #))+ А+ (х, :) + отт [ем (х,2кл де : (12) 


справедливых для О х<А,Ё>0, где Г = {2 : -№ + ,, ю+ 0%, | — контур интегрирования в комплекс- 
ной плоскости, представляющий прямую линию, параллельную мнимой оси и отстоящую от неё на ве- 
личину а, в которой &4 выбирается таким образом, чтобы все полюса подынтегральных функций в 


(10)-(12) были бы левее @Ё,, 


К(7)= ( _ у?) 2 5Ну2-у Ин 2 спур — у $Н/Р $Ну2 — 1), (13) 
№ (х,2) = М2 -у?2)н у2 $"(/2хи-), (14) 
М№М,(х,2) = ( Е у?2)св Утзнуг у Инт хи св у2-— у\/2 $Н Утэнгжи“)- сп(у2(^ _ х)в"), (15) 
М. (х,2)=У 2№ (х,2)5Пу2 - сн(угхй (2) (16) 

№(х,2) = ( = у? }н/ рану + у Иу 2 6н/2енугм-")- $н(/тхи-'}нуг)- $®(у2(^ = х)в"), 
К, (2)= $1у28(2), (17) 
ЛиБ(2) = 2 Ц2ехр(- 26 „)-1)+ (2-Е) '(- ехр(- (2 -=6 к), (18) 
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об -Е- Ах, 





где а -— постоянная из (7), а -из (8), безразмерные параметры у и и определяются формулами 
ук у. 6 жил, 
ай К 1-2и 
Внеинтегральные слагаемые получены в формуле (12) в результате выделения обобщенных квадратур 
[14-16] с помощью выражения Р(2)сН(у2хй`!)/5Нух . 
Квадратуры в формулах (10)-—(12) существуют при выполнении условия Ё> (р-х)/а и благо- 


даря алгебраическому убыванию подынтегральных функций на бесконечности 


В-х 
Б(2)М:(х,2)Е-К2) = ие при |2 |-> ®°, (19) 
в г В-х м В-х 
р [се о й а я] 1 й 25" . 
М (х,2)Кз (2) 7 ры (20) 


при |2 |-> ®. 

Исследование подынтегральных функций в формулах (10)-(12) показывает, что все они меро- 
морфны в комплексной плоскости переменной интегрирования г = & + Я, то есть имеют в качестве 
изолированных особых точек только полюсы, которые доставляются обращением в ноль знаменате- 
лей этих функций: А(2), К,(27). Кроме того, следует отметить, что подынтегральные функции в 


(10)-(12) при 2=0 и2=уУ”? обращаются в ноль. 


№-(х,0) = М, (х,0) = М(х,0) = ©) =, (©)=0, 
(21) 
М (ж, у?) = Мож, у) = МС, у?) = у) = у )= 0, 
причём Л. (х,2) и К (2) имеют двукратный ноль при 2 = 0. 
В заключении этого пункта заметим, что выделение главной части поведения подынтегральной 
функции в формуле (12) для о(х,Е) производилось с помощью трансформанты Лапласа решения со- 


ответствующей одномерной упругой задачи о внедрении недеформируемой полуплоскости Вв упругую 
полосу на жестком основании. 
Нули функции ^(7) в комплексной плоскости. Для вычисления контурных интегралов в формулах 


(10)-(12) методами теории функций комплексного переменного необходимо знание нулей функции 
А(2) и их свойств в комплексной плоскости. Нули функции А(2) из (13) определяются в комплексной 


плоскости 2 =& + из решения трансцендентного уравнения 
К(7)= ( _ у? 2]в 2 зпу2 -у Ис \2 спу2 — УР $Н/Р $Ну2 — 1)= 0, (22) 
которое совпадает с соответствующим характеристическим уравнением [5]. 

При численном определении корней уравнения (22) в комплексной плоскости 2 = & + М исполь- 
зовались итерационные численные методы определения корней, требующие хорошего начального при- 
ближения. Для получения начального приближения корней уравнения применяется параметрический 
анализ функции (2) по двум безразмерным параметрам у и и. Зафиксировав у и положив И=Ов 
(22), получим упрощенное уравнение 

( зы у? н/Ну2 =0, (23) 
решение которого распадается на решение трех уравнений: 
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1-\У27=0, 61/2 =0, $Пу2 = 0. (24) 
Первое из уравнений (24) дает корень 2 = У”, являющийся корнем не только уравнения (22) 


для любых у и И, но и всех остальных подынтегральных функций, что было отмечено ещё в (21). 


Из последнего следует, что 2 = \`* не является полюсом подынтегральных функций в формулах (10)- 
(12). Второе и третье уравнения в (24) приносят две бесконечные последовательности нулей: 


2 
П 
Рок = 2 (1 в 2) 7! К = 0,1,2, = #1, (25) 
2 = = ИЗ (26) 


Нули 2х (25) и 2%. (26), определенные при И=0 из уравнения (22), являются начальными 
приближениями для последовательного с увеличением и определения соответствующих по номеру 
нулей 2% = И), ЕО. М, = 281), п = 1,2... из численного решения уравнения (22) при 
фиксированных у. 

Нули 2^ =2,(0) (25), соответствующие И= 0, располагаются на отрицательной части дей- 
ствительной оси Т1т(2) =п =0 комплексной плоскости 2 =&+ я. Множество 2 при И = [0,®) пред- 
ставляет полубесконечный промежуток [20(0), ©) действительной оси, где 2%(0)= —п2/4, апри И=2 
и произвольном у 2о(И) попадает в начало координат 2,(2)=0. Таким образом, при 0<И<2 
выполняется неравенство =? 14 2о(и) <0, а при И>2 значения 2о(и) >0 и располагаются на 
положительной части действительной оси (рис. 2). При этом 2о(и) единственный из всех корней (22), 


который при произвольном ‘у обладает свойством т 2,(И) =®. Множества 2, ), К =1,2,3,.., 
И>о 


И = [0, ©) при произвольном у представляют собой отрезки на отрицательной части действительной 
оси, левый край которых совпадает с 2,„ = 2,(0)<0, К=1,2,3,..., из выражения (25), а правый пред- 


ставляет точку сгущения йп 2%(И) = 2, (®)> 2, (0), где 2% (0) < 0 (на рис. 2 стрелкой указано направле- 
И>о 


ние расположения точек множества 2» и) при И > ©, совпадающее с положительным направлением 
действительной оси). 

Заметим, что при больших значениях А 2% = 2к(0) из (25) имеют очевидную асимптотику 
2к(0)= —(п^)} +0) при А->ю. 


Нули 20) =2#(0) (26), соответствующие И =0, располагаются на мнимой оси Ве(2)===0 
комплексной плоскости 7 =&+ Я, являясь начальными приближениями (при фиксированном у) для 


^ 


последовательного с увеличением И определения в комплексной плоскости соответствующих по но- 
меру нулей 2# и) из численного решения уравнения (22). Ноль 2% = 21(0) = 0 не является полюсом с 


учетом соотношений (21). При у=1 множества 2:0), И Е [0, ©) представляют в комплексной плос- 
кости г =&+ м траектории в виде кривых линий, похожих на незаконченные эллипсы (см. рис. 2) с 
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главными осями, параллельными действительной оси 1т(2)=п=0. С ростом л эллипсы располага- 
ются периодически с периодом п (у=1) выше и ниже действительной оси, причём их ширина с ростом 


номера уменьшается (рис. 2). Уменьшается с ростом номера /7 и ве(== (7) = ==, то есть для корней 
2:7) уравнения (22): 

т Ве(2* (= т Е =0. (27) 
—>0 По 


п 


^ 


На графиках 2:7) стрелками указано направление расположения нулей с ростом аргумента И ‚ при- 
чем при возрастании И от 0 до И,, действительная часть 2) возрастает и при И =И., достигает 


своего наибольшего значения. При дальнейшем увеличении И от И,, ДО <х величина ве(2= 7) 


убывает до некоторого конечного значения Ве(2* (х)} (рис. 2). 





















































Рис. 2. Расположение нулей А(2) в комплексной плоскости 2 при И е [0,°) , у=1. 


В левом верхнем углу — общий план расположения нулей, в нижнем левом углу — 
расположение множества точек 2(И) в увеличенном формате 
м я 2 22) ие = п: 6 ыы 
нимая часть корней 2»; уравнения (22) 1т\2;„И |= п; с ростом л приближается к + г ? 
что также следует из уравнения (26): 
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ии(2 (7 = ПЕ = Ра: +0(1), п=1»2,,... (28) 
В этом случае для больших номеров /7 нули 27) можно записать в виде 
ры ити, и (29) 


где Ш" & =0 из (27). Подставив (29) в уравнение (22) и считая при этом & малой величиной, для 
П—с 


больших л получим следующую оценку: 


в - +6[ при Л —> <, (30) 
у уп п 


а корни 2:7) при этом для больших значений / получат асимптотику 


р -<[ +120) при 7->®. (31) 
ууу (7 


Следует заметить, что графики 2/ и 2, обладают симметрией относительно действительной 








оси 1т(2) =п=0, как и подынтегральные функции в квадратурах (10)-(12), а для 2х выполняется 
равенство 

ри (32) 
На всех траекториях нулей (и) п =1,2,3,..., Указаны точками значения 2,(И,) и 2 где 
= 1,2,3, причём и =1, и, =2, и = 2,7. Значения 2), п=1,2,3,..., при одинаковых и, ‚= 1,2,3 
соединены пунктирными линиями. 

Графики нулей функции К2: 2-2), П.К = 1,2,3,..., 2), представленных на рис. 2, 
рассчитаны с помощью программы, написанной для системы Маре, при существенно увеличенном ко- 
личестве значащих цифр в расчетах. 

Анализ полученных решений. После вычисления контурных квадратур в комплексной плоскости в 
(10)-(12) [17, 18] для функций 7(х,Ё), и(х,6) и о(х,Р) получаются удобные для вычисления фор- 
мулы в виде рядов по полюсам подынтегральных функций: 


Т(х, О Их ВО], аи. 250. (33) 
ик) = во 5х, |, 0<х<А, Ё>0, (34) 
К=1 
обв) = МИ № ов (6 з(х, 6+ к, (35) 
с ма в. ера. №5; 
К=1 
где $„(х,Ё) (К= 1,2) выражаются формулами: 
$ (х,6)=5(х, 6) +б(х,Ё), (36) 
С(х,Ё)=28е Уве, х, Е) +56 ЖЕ), (37) 
К=1 К=0 
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$, (х,Ё) = $ (х,Ё)+б(х,Ё)-26(х,Ё-Ёк), (38) 
5 (х, Е) = КЕ, х,Ё)-К(0,х,0), —55(х,Ё)=К(0, Хх), 
К(2,х,Ё)= м ехр(2Ё), К“’2,х,Ё)= м ехр(2е), 
(хо) Пт”, 62) ЕЕ 
5°(х,Ё)= ое об дих, (39) 
За 5х, #-)-50(х,Ё), (40) 


где А’_(2) есть производная по 2, //(Ё) — функция Хэвисайда. 
В формуле (35) под А(Р) понимается нормированная функция Д(Ё)= =е" (Е, - ВН(Е) из (4), 
где &=и,/Д,. Функция б(х,Ё) (37) в формулах для $,(х,Ё) (К=12) содержит два бесконечных 


ряда, один из которых по полюсам 2%, К =1,2,3,..., следующего вида: 





5*- Ук, Ых,2)-0 Мх,2). 41 
Узор, ве) (4 


сходимость которого не очевидна, так как Ке(2,)>0, Ё>0. С учетом асимптотики 2х (31) при 


и 
РЁ 


К —> ® можно показать, что последовательность экспонент {в | является ограниченной ДлЯ фик- 


сированных значений # ©- ТЕ =]: 


ей |< М (м>0), 








а ряд (41) является сходящимся [19], в том числе за счет убывания 6(2) при || > ©. С другой 


стороны, учитывая упорядоченность Ве(2: )}. где для двух соседних членов выполняется условие 


Ве(2; }> Ве(2й.), К =1,2,3,..., 
сумму ряда (44) можно представить в виде 
2 [5+ 7+7” +7 Ь г 
я в || + Ув ( 2} ры | В(х,2)= д (42) 


где ряд в скобках является сходящимся, причем со скоростью геометрической прогрессии, так как 


Ве(2; ыы 0, А=2,3,... Из выражения (42) вытекает, что при неограниченном возрастании Е 

(Е —> оо) , неограниченно возрастает 5” из (41), а решение, его содержащее, является неустойчи- 
вым при любых значениях И>0. 

В формуле для С(х,Ё) (37) из $,(х,Ё) (К=12) , кроме вышеупомянутого ряда (41), содер- 

жится и другой бесконечный ряд по полюсам РКО: Учитывая свойства 

2, изученные выше, а именно то, что при 0 < И <2 полюс дИ)< 0, при И =2 полюс 2@)= 0, при 


И>2 полюс д И)> 0, можно утверждать, что ряд по полюсам 2, 
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5 = Ур, : (43) 


К=0 
где Ь(х,2) из выражения (41) для фиксированных значений Ё сходится со скоростью геометрической 
прогрессии для любых у и И. Однако при И>2 ‚ когда 2> 0, сумма ряда 5” (43), перестро- 
енная по формуле 


5 = о Н Ух, Не : (44) 


К=1 
где В(х,2) из (41), при Ё —> ® неограниченно возрастает, хотя ряд, стоящий в скобках, сходится. 


Это означает, что при И> 2 решение, содержащее ряд (43), является неустойчивым. Таким образом, 
решения рассмотренной задачи (33)-(35), содержащие ряды (41) и (43), являются термоупругонеустой- 
чивыми, начиная со сколь угодно малой скорости скольжения И > 0. 

Численный анализ полученных решений. Для детального изучения термоупругонеустойчивого ре- 
шения рассматриваемой задачи используются формулы (33)-—(35). Числовые значения параметров за- 
дачи следующие: [м = 44 ' 109 Н/м>2, у = 0,25, а = 1,04 : 102 КТа=4`: 103 м/с к=8`: 104 м?/с, 
К = 62,8 Вт"м/К (параметры материала соответствуют табличным значениям для чугуна), / = 0,15, 
Ло =олА, и =0,01 м/с, А=2-103 м, И= 1,58 м/с, & = 6,93-103с, Ё, =5.1077 с, в результате чего 


безразмерные параметры задачи у и и приобретают следующие значения: \у =10.10“ и и = 0,25. 
На рис. 3 представлены графики изменения смещений 1/(х,#) по толщине покрытия (0 <х< В), рас- 
считанные по формуле (34) при разных значениях х, которые обозначены цифрами 1—5, причём 1 
соответствует х= 0, 2-0,25/, 3-0,5/, 4- 0,754, 5 - А (на контакте). На рис. 3 видно, как со 
скоростью геометрической прогрессии развивается амплитуда смещений неустойчивого решения 
и(х,Е) на собственной частоте 1т(2!), теряя физический смысл с некоторого момента времени, 
причём при х = 0,5А это происходит немного раньше, чем при других х. 

На рис. 4 приведены графики изменения напряжений о(х,Ё) для различных глубин по х, кото- 
рые, как и на рис. 3, обозначены цифрами 1-5. Начиная с некоторого момента времени, амплитуда 
колебаний о(х, г) на собственной частоте 1т(27!) нарастает по экспоненте и неустойчивое решение 
теряет физический смысл, причем примерно в тот же временной период, когда это же происходит со 
смещениями и(х,Ё) на рис. 3. С помощью врезки на рис. 4 подробно показано изменение величины 
напряжений о(х,гГ) на отрезке Ё е› [0, 1,8.10°] в более мелком временном масштабе: различается 
запаздывание по времени прихода волны сжатия в сечение х, наибольшее при х = 0 (отмечено цифрой 
1), амплитуда пришедшей волны сжатия при х = 0 удваивается за счет отражения от жесткого основа- 
ния (две цифры 1 на одной и той же прямой), двумя цифрами 5 отмечается рост амплитуды напряжений 
сжатия на контакте при х = Аит.д. 

На рис. 5, а представлены графики температуры 7(х,г) в покрытии для различных значений 
х, обозначенные цифрами 1-5, как на рис. 3, 4. С помощью врезок 6, в, гдаются наиболее существенные 
фрагменты графиков температуры 7(х,Ё) в более мелком временном масштабе: на врезке 6 даны 
графики развития температуры на начальном временном отрезке при различных значениях х, график 
с цифрой 1 соответствует выполнению при х = 0 граничного условия 7 (0, 2) = 0, с увеличением хгра- 
фики 7(х,Е) показывают всё более быстрое нарастание температуры с максимумом на контакте при 
х= А (кривая 5); на врезке в представлен график температуры 7(р,Ё) (на контакте) на временном 
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отрезке Ё > 3. ‚ показывающий колебательный характер графика на собственной частоте 1т(2!); на 


врезке гпоказан интервал времени, на котором амплитуда колебаний претерпела 50-кратное увеличе- 
ние по сравнению с врезкой В. 
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Рис. 3. Графики распределения смещений и(хд) по толщине Рис. 4. Графики распределения напряжений с(х,{) по толщине 
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Рис. 5. Графики температуры 7(х,д на разных глубинах по толщине 
покрытия и на различных временных интервалах 


В дальнейшем (рис. 6) можно видеть экспоненциальное развитие амплитуды колебаний графика 
КА,6), достигающей десятков градусов, что свидетельствует о потере физического смысла для х = А 
при & > 0,18 с, для х =0,75А- при & > 0,82с, для х=0,5А —при Ё > 1,42 с, для х =0,25й- при Ё > 2,05 с. 
Такая зависимость нарушения физического смысла оказалась свойственной только температуре, опре- 
деляемой из параболического дифференциального уравнения. 
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Рис. 6. Температура Х(х,ё на разных глубинах по толщине покрытия с указа- 
нием времени потери физического смысла 


Выводы. Термоупругодинамическая неустойчивость скользящего фрикционного контакта при движе- 
нии жесткой полуплоскости по упругому покрытию на жестком основании наступает при любой скоро- 
сти движения полуплоскости. Неустойчивость решения рассмотренной задачи о скользящем фрикцион- 
ном контакте заключается в том, что, начиная с некоторого момента времени, полученное решение 
теряет физический смысл — амплитуда собственных колебаний неограниченно возрастает по времени. 
Следует заметить, что решение задачи для смещений и напряжений теряет физический смысл раньше, 
чем решение для температуры, что может быть существенным моментом при разработке методов диа- 
гностики возникновения термоупругодинамической неустойчивости скользящего контакта. 
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\. В. РаепЕ5ом, В. Г. МИит, $. $. МоКом, А. $. УазПуеми 


А опе-атепзюпа! пептоейв$йс сотасЕ ргоМет оп те иегйса/ п5егйоп оЁа По ВеЕ-р/апе тоутд Вопгота/у аЁ а 
сопапЕ зрееа оуег те ейа5йс соайпд (тр) ие те БоНот $/ае о! {те а@ег по@\У гезёпд оп пе поп-деюттта 
Юипаайоп 15 сопяЕГед. Оп {те Гоипаайоп зигтасе, Пе {етрегайиге 15 КерЁ сопЯапЕ. А ВеаЁ Пои/ депегаеа Бу пе 
свопа/ сопвасЕ [5 Чгесеа ю пе соайпо. Тре ргоет зо/ивоп (5 оатед итд те [ар/асе Иедга! гапзЮюпт апа /5 
Гергезеже4 т {те югт оЕсотоиг тедга(5. Тре юсавоп о {те зоивоп тедгапа ро/ез 15 5/4 аЁ иатоих ва5К орйопе. 
Тетрегайиге, @зр/асетепЕ апа 525$ @5трийоп$ оуег те соайпд аерЁ? аге депуеа т {1е отт о! {те пПИпЁе епез 
оуег едеп'ипсвопе. 1: 15 звоип таЁ те тегтоеазюаупатис 15 ау о! {пе оматпеа зо/инопз 15 ргезепЁ асго55 пе 
ипо!е Ите иегиа! апа аЁ апу уе/осйу о {те ра/!!-р/апе $/пд оуег йе соайпд зигТасе. 

Кеуигога$: тептое/взюаупатис пару, соцр/еа {пегтоея$йсйу рго ет, тсвопа!/ сотасЕ, соайпд. 





* Тне гесеагсн {5 Чопе мн Не Япапсиа! зиррог* гот ВЕЕГ дгапз пос. 14-08-91166-СРЕМ_а, 14-07-00271-а, апа 14-07-00705-а. 
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